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OS : 

Q\ , We suggest treating a conducting network of oriented polymer chains as an anisotropic fractal 

whose dimensionality D — 1 + e is close to one. Percolation on such a fractal is studied within the 
^q' real space renormalization group of Migdal and Kadanoff. We find that the threshold value and ah 

r—{ , the critical exponents are strongly nonanalytic functions of e as e — > 0, e.g., the critical exponent of 

conductivity is e~ 2 exp( — 1 — 1/e). The distribution function for conductivity of finite samples at the 
percolation threshold is established. It is shown that the central body of the distribution is given by 
a universal scaling function and only the low-conductivity tail of distribution remains e-dependent. 
Variable range hopping conductivity in the polymer network is studied: both DC conductivity and 
AC conductivity in the multiple hopping regime are found to obey a quasi-ld Mott law. The present 
(— i ■ results are consistent with electrical properties of poorly conducting polymers. 

9 ' 

• $H ■ I. INTRODUCTION. 

: 

Charge transport in structures with fractional dimensionality has attracted a high degree of attention due to both 
its fundamental and applied interest |l],^] . Conductivity of a random fractal of resistors was studied in Ref . || as an 
example of critical phenomena. A large body of literature has been devoted to a study of the unusual dynamics of 
electrons in regular and random fractals . The discussion of the questions and appropriate references may be found 
£h ' in recent reviews (see, e.g. Jill). 

In the present paper we consider problems of percolation and hopping transport on nearly one- dimensional strongly 
anisotropic fractals [Q with dimensionality D = 1 + e. These fractals are expected to exhibit unique properties because 
their dimensionality at e -C 1 is close to the low marginal one for the percolation transition. In contrast to isotropic 
fractals a small parameter e enables us to obtain the exact solution. Moreover it becomes possible in the case of 
nearly-one-dimensional fractals to establish not only average characteristics but their entire distribution functions. 

Another motivation for a study of quasi-one-dimensional fractals is recent experimental data on conducting polymers 
such as doped polyacetylene, polypyrrole and polyaniline [^|-^2|. In general, this class of polymers has a great variety 
of transport properties. In the undoped state these polymers are semiconductors with an energy gap of Peierls-Mott 
■ origin |U5). With doping the energy gap is suppressed quickly and for the highly doped case there is finite density 
of states at the Fermi level. The room temperature conductivity (ctrt) of heavily doped sample may attain metallic 
ON ■ values and temperature and frequency dependencies of conductivity may be close to metallic. The nature of the 
metallic phase in these samples is presently a subject of intensive study. One point is that the metallic state is 
dependent upon by strong interchain coupling |Io|- |r^ , |l4f . 

In doped polymers with moderate grt (of the order of several hundreds S/cm) the conductivity, as a rule, decreases 
with decreasing temperature ||,[L(^-[L2|. Because this decay follows a power law in a large temperature interval, 
presumably, these materials are near the metal-insulator transition which happens at the critical interchain coupling. 
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Poorly conducting doped samples with ort of the order or less 1 S/cm have behavior that can be classified as 



"dielectric" 10 : it is similar to that observed in amorphous semiconductors. For such materials, DC conductivity is 
strongly dependent on temperature and, generally, follows: o~dc ex P — (Tq/T) 1 / 2 . It is noted that caution must be 
• i-H , made to review a large temperature range in comparison of experimental conductivity with the model dependencies 

X: & 

For a variable range hopping (VRH) mechanism of transport the temperature dependence of conductivity was 



initially derived fll5| , [l6| to be : udc °^ exp — (To/T) 1 ^ d+1 \ where d is system's dimensionality. For d = 1 this formula 
reproduces the observable dependence o~dc °t exp — (Tq/T) 1 / 2 . However this approach is not correct, since the Id 
VRH |nj yields the Arrhenius law, croc oc exp —Tq/T with Tq set by the highest barrier that occurs in the chain. 
The Id VRH was modified for a quasi-ld system (TJJl^] to include weak hops between the nearest-neighbor, thereby 
avoiding the highest barriers. This approximation results in a quasi-ld VRH law, aoc oc ex P — (Tq/T) 1 ^ 2 . 

Experimental measurements of microwave conductivity and dielectric constant in poorly conducting doped samples 
[ p~0| , p~T] | revealed that both are strongly dependent upon temperature too, most probably according to the same quasi- 
ld Mott's law, i.e. exp — (Xq/T) 1 / 2 . The usual theory of hopping transport predicts, however, only a very weak 
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power temperature dependence for the frequency-dependent conductivity and the dielectric constant in two- and 
three-dimensional systems jl9| . 

In the present work we exploit the specific structure of the polymer network to understand these peculiar features of 
conducting polymers. In stretched polyacetylene this network is formed by coupled polymer chains oriented along some 
direction. Electron micrographs shows that in these substances polymeric chains are organized into fibrils , which 
may be distinctly seen to be subdivided into smaller ones |2Q| . In a non- fibrillar form of conducting polymers, like 
polyaniline, X-rays data reveal the existence of highly ordered "crystalline regions" with metallic properties [p~p|— p~2|| . 
Therefore the whole network of stretched polyaniline may be thought of as constructed from long one-dimensional 
polymer chains randomly coupled by metallic islands of various sizes. The volume fraction of metallic islands can be 
small. 

We assume here that polymer structure represents a nearly one- dimensional fractal. That means a specific kind 
of polymer chain organization, defined in the following way: Choose a three-dimensional cube with the edge L. 
Chains, which are coupled within this cube, form a set of bundles disconnected from each other. If for large enough 
L the cross-section of the maximum bundle is proportional to L e , where < e < 2, then we shall call the system 
d* = 1 + e -dimensional. Obviously e = for purely one-dimensional systems (sets of uncoupled chains). Note, 
that if one assume chains to be connected either with a low concentration of uncorrelated interchain links, or with 
weak links (their resistivities being high compared to intrachain ones in our example), then we are dealing with a 
quasi-one-dimensional system ]l4fl , which is three-dimensional according to our definition. 

The problem of electron localization in similar fractals was studied in Refs. |2l],|2^]. It was found that even in the 
presence of a weak disorder all the electronic states remain localized as long as e < 1. Therefore, a mechanism of 
charge transport in the fractal with e < 1 is supposed to be variable range hopping (VRH). This assumption is in 
agreement with the experimental observations for poorly conducting highly doped polymers for which there is a finite 
density of states at the Fermi level |2j| . 

The usual method to treat VRH models is the effective medium approximation p9| , which gives wrong results in 
the nearly-ld case. For example, for the percolation model this method gives the threshold concentration of broken 
bonds ct ~ e, while ct ~ exp (— 1/e), as we shall see later. The results for critical exponents are also wrong in this 
approximation. To treat VRH in a nearly one-dimensional fractal we choose the following approach. We will first 
study the percolation problem in a nearly-ld fractal exactly. The VRH model is reduced to the percolation problem 
by constructing the effective percolation lattice [ p^|j2^j25| ] . 

In this way we have found that at low temperatures the VRH conductivity obeys a quasi-ld Mott law: (Jdc °^ 
exp — (Ti/T) 1 / 2 but the characteristic temperature T\ is greater than To for Id chain by a factor 1/e. Similar temper- 
ature dependence is obtained for AC conductivity. These results can explain the observed temperature dependence 
of conductivity and dielectric constant in poorly conducting polymers. Additionally it was shown that there is the 
strong frequency dependence of conductivity in the region of extremely low frequencies. These peculiarities reflect the 
fact that in the random fractal with dimensionality close to one the low frequency conductivity is entirely controlled 
by the weak charge transfer between clusters. Each cluster is very dense and remains well isolated. 

There exist several different problems related to the percolation. First, one can be interested in statistical properties 
of percolating media: distribution of connected clusters, probability of two or more points to be connected, etc. As it 
was shown by Fortuin and Kastelein this problem can be reduced to the g-component Potts model in the limit 
q — > 1. Thus, the powerful set of field theory methods may be applied. This analogy, however, does not allow us to 
treat the conductivity of a percolating cluster. The evaluation of the conductivity exponent /i, which describes the 
DC conductivity a behavior near percolation threshold: 

is much more complicated task than the "field-theoretic" ones, such as the exponents of correlation length and of 
infinite cluster capacity, etc. 

Thus our first aim is to study the critical behavior of conductivity near the percolation threshold in a d-dimensional 
lattice, where d is close to lower critical dimensionality, i.e. d = 1 + e, e <C 1. The real space renormalization group of 
Migdal and Kadanoff (RGMK) p7[-pQ], being exact at d = 1, may be expected to be the appropriate tool as d tends to 
unity. This method was applied to the percolation problem several years ago by Kirkpatrick |3|] . He had found critical 
exponents of correlation length and conductivity by using the RG equations for conductivity distribution truncated 
up to the first moment. Though he had not considered explicitly a case of nearly-ld system, this method, if properly 
applied, gives the right dependence of conductivity exponent on e except for a pre-exponential factor. 

We extend the RGMK method to consider the conductivities and resistivities distribution functions in random 
media. This enable us to derive the equation for the conductivity exponent /i, which gives realistic values not only 
in the nearly-ld case. Moreover, explicit expressions for the distribution functions at the percolation threshold will 
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be obtained. To the best of our knowledge only numerical estimates of random conductivity momenta were available 
until this work (see, e.g. |3l|] ). 

Thorough investigation of the percolation problem and its various modifications is of interest both for its conceptual 
significance and due to its numerous applications |]6|, j32||3"3f| . Beside its application to random conducting media on 
which we shall concentrate here, the percolation approach was used, e.g., to treat the rigidity transition in random 
networks [j34| , and mechanical breakdown in solids [p5| . Another application is the magnetic flux flow in type II 
superconductors fl36[ . In the case of magnetic vortices pinned by disorder, their motion just above the depinning 
threshold happens along a sparse (possibly, fractal) network of persistent channels. Directed percolation model is 
often applied now to describe a wide class of phenomena, in particular, self-organized criticality |37|j3^ ]. 

The paper is organized as follows: In Chapter [n] the notion of fractional dimensionality is introduced for the oriented 
chain arrays and illustrated by hierarchical structures. In Chapter III RGMK transformations for conductivity of 
disordered media and for connectivity of percolation system are derived, the latter one is studied in Chapter IV. 
The RG equation for the distribution function of conductivity at the percolation threshold is solved in Chapter [V, 



and the explicit form of distribution function is found in Chapter VI. Using sca ling relations, results for the AC 
conductivity of the percolating lattice near threshold are obtained in Chapter VII , and they are applied to describe 



the temperature and frequency dependence of conductivity for the variable range hopping transport in Chapter VIII 
Results are discussed in Chapter IX. Three Appendixes contain technical details. 



II. HIERARCHICAL COUPLING OF ORIENTED CHAINS 



The (m, n) hierarchical structure (HS) is constructed through the infinite repetition of two successive steps (see 
Fig. [j]): a) construction of n-chains, and b) construction of m-bundles. After every step the resulting construction 
may be treated as a new bond (l-th level bundles). 

We use the following definition of dimensionality for a chain fractal considering an array of one-dimensional chains, 
connected in some regular or random fashion by transverse bonds of various lengths. In an L-size cube chains form 
a set of bundles, connected inside this cube. Within each bundle in the cube the chains are interconnected. There 
are no connection between the bundles within the L-size cube. If the number of chains in the maximum sized bundle 
scales as L £ for large enough L, where < e < 2, then we have D = 1 + e -dimensional network. Obviously e = for a 
purely one-dimensional systems (sets of disconnected chains). The characteristic feature of the fractals, constructed 
from oriented Id chains is their self-similarity: the system at any scale looks like subdivided into bundles, which in 
turn are subdivided into smaller ones, etc. 

In particular, the dimensionality D — 1 + lnm/lnn may be ascribed to the (to, n) hierarchical structure, if in spatial 
dimension d > D we replace every site with 2m l bonds attached to it (I is the level of bundles attached at each side 
of this site) with m l sites connected by transverse bonds of infinite strength (see example in Fig. |^). 

Our hypothesis here is that oriented polymer network structures are of this type (with D = 1 + e close to 1, e <C 1), 
at least in some wide enough interval of length scales, e.g. from the scale of fibrils (hundreds of nm) down to molecular 



scales. Transmission electron micrographs of fibrillar polyacetylene (see e.g. 1 20 ) appear to support this hypothesis. 



Of course, the real structures are not regular ones, and the requirement of self-similarity here is to be treated in 
statistical sense. Nevertheless, we shall use the RGMK scheme, based on regular fractals (HS) for their analysis. In 
the case of conducting polymers such as nonfibrillar doped polyaniline and doped polypyrrole one may assume the 
polymer networks apparent fractality in some scale range to be caused by a dilute distribution of crystalline regions 
providing interchain links (fractality, generated by randomness |39[| ). Having in mind that the RGMK is exact in one 
dimension, one may hope to obtain meaningful results when the dimensionality is close to 1. 



III. MIGDAL AND KADANOFF EQUATIONS 



The renormalization group of Migdal and Kadanoff may be formulated in a quite simple phenomenological fashion. 
Suppose we have some random D -dimensional medium with fluctuating local conductivity/resistivity. Let us consider 
a A-size cube within the medium. Its conductance is a(\)\ D ~~ 2 , where the random conductivity a (A) is A-depcndcnt 
for strongly inhomogeneous systems. The distribution function for conductivity and for resistivity p(A) = l/er(A) are 
defined in the Laplace representation as: 

P(ry, A) =(exp (-r,a (A))); Q(s, A) = (exp (-sp (A))) . (2) 
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If we change the size of the cube, A — + A' = nA, we arrive at some new random variables cr(A'), p(X') with distribution 
functions P (n, A'), Q(s,A'), respectively. The cube's enhancement may be treated as n-times expansion in one 
( "longitudal" ) spatial direction, and in D — 1 other ("transverse") ones. If one intends to treat these transformations 
as infinitesimal ones afterwards, the order of operations is not important. 

The RGMK scheme is based upon two approximations: i) enhancing the size n times in the longitudinal direction 
is treated as connection of n resistors in series, and ii) in a similar way, the transverse cube's enhancement is replaced 
by the parallel connection of m = n ^ 1 elements. Thus we have: 



n 1 m 

(A) = -$>(A) , a(nX) = -J2* { l n) W ■ 

11 — ' m. — ' 



(3) 



Here the tilde values refer to the rectangular element with the dimensions nA in the longitudinal direction and A 
in other ones. In both steps the resistivities pi and conductivities 07 = l/j5j of the constituent component are 
supposed to be independent random variables, and, therefore, we have the distribution functions to be transformed 
in two steps simply as: 



Q(s, A) = (exp (-spW (A)) ) - Q n (s/n, A) , P (77, nA) = P m (n/m, A) , m = n 



D-l 



(4) 



This transformation is exact for the (m, n) hierarchical structure. 

Equation (Q) should be supplemented with the relation between conductivities and resistivities distribution functions 
(DF) in the Laplace representation. It may be easily derived from the definitions (||), using the following integral 
identity: 



-x/a 



1- Vi / 
Jo 



dy 



where J\ is the Bessel's function. As a result, the relation between Q(s, A) and P(n, A) takes the form of Hankcl's 
transformation: 



Q(s,X) = l 



dr) 



si — J 1 (2 v ^)P(r?,A) 
/o Vv 



(5) 



the reverse relation is of the same form. This transformation has the following properties to be used later: 

Q(0,A) = l-P(+oo,A) , Q(+oo,A) = l-P(0,A) , (6) 



dQ(s,X) 



dr] . . dP(i],X) 

lo VV dr > 



(7) 



d 2 Q{s,X) 
ds 2 



•s / ^J 1 (2Vslj)r 1 P(7 1 ,X) 
h Vv 



(8) 



Thus we have a closed set of equations for arbitrary rescaling factors n and m. It appears to be more convenient to 
deal with an infinitesimal transformation by setting n = 1 + <5A/A, m = 1 + eSX/ X. From Eq.(^) we have the variation 
of distribution functions consisting of longitudinal and transverse parts: SP = <5;P + S t P, SQ = 5iQ + 8 t Q. 



s dQ(s,X) + Q (S)A)lnQ(SiA) 



ds 

d P(r,,\) 
drj 



P(n, A) In P(n, A) 



SiQ(s,X) = 
5 t P(i], A) = e -11- 
Using the relations (^J?]), one may rewrite the first of these equations as: 
SiP(vA) 



sx 
T ' 
sx 

A ' 



v dP(ri,X) _ r^ Jl(2 ^)Q( S) A)lnQ( S ,A) 
°V Jo V s 



SX 

T 



(9) 



(10) 



Then, from SP — §x<5A = 5[P + StP, we have the following equation: 
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(1 - e) V dP %' A) + eP(n, A) InPfa A) - f° Ji (2^) Q(«, A) lnQ( S , A) . (11) 

This equation, combined with the Eq. (|^), determines the evolution of the distribution function upon size scaling in 
a closed form. 

This scheme also allows us to treat the percolation system, introducing the probability that the A-sized cube is 
disconnected (i.e., has zero conductivity or infinite resistivity) c(A). Taking into account the definitions of distribution 
functions (^|), c(A) may be written as: 

c(A) = P (+0, A) = 1 - Q (+oo, A) . (12) 
Putting in Eq. ( |Tl| ) = +oo and using the formula (^), we have: 

dc 

— = eclnc-(l-c)ln(l-c) . (13) 
gA 

The right hand side of this equation has three fixed points: two stable ones, c — and c = 1, corresponding to 
connected and disconnected systems, respectively, and the unstable fixed point, c = q, < c t < 1, 

ec t \nc t = (1 - ct) In (l-ct) , (14) 

corresponding to the percolation threshold. 

Now let us consider statistical properties of clusters for the percolation problem, i.e., the distribution of clusters 
over sizes and site numbers, existence and capacity of the infinite cluster, etc. We suppose every bond of a HS, to be 
either broken, with probability c, or not, — with probability 1 — c. The statistics of percolating network is closely 
related to the thermodynamics of the g-states Potts model p6f . Namely, if we consider the partition function of the 
latter: 

Z q = exp {-H q ^ , (15) 

hW=kJ2(i-Kvi), (is) 

(ij) 

where variables r)i = 0, 1, . . . , q — 1, 8 is the Kronecker (5-symbol. Z q may be expressed in terms of the percolation 
model as: 

Z q = (q T ) ■ (17) 

Here T is the total number of connected clusters in the percolation model, and the average is over realizations with 
the broken bonds concentration c = exp(—K). 

To establish a further relationship, it is necessary to introduce external fields into the Hamiltonian of the Potts 
model: 

<'=«f+''iE( 1 -U' ( 18 ) 

i 

which may be thought of as additional bonds with strength h\ between every site of a given lattice and some fictitious 
external one. 

The statistical properties of the percolation model may be obtained from its "free energy" 



f(K,h 1 ) = ~-^-lnZ g (K,h 1 : 



= -^< r >/u, (19) 

9=1 



where N is the total number of sites. For an example, the order parameter 



P = l- |f , (20) 

1 hi=0 



G 



characterizes the "capacity" of the infinite cluster, i.e., the probability of a given site to belong to the infinite cluster. 
The second derivative is the "susceptibility" : 



dh\ 



(21) 



hi=0 



which gives the average number of sites in finite clusters || . 

To obtain the RG equation for / of a HS, one should sum over intermediate sites in bundles, thus, performing a 
transition from the initial Hamiltonian, containing variables of all sites, to the one, containing variables of sites of 
the next level (see Fig. ||). This transformation reproduces the structure of the initial Hamiltonian ( |l8| ) with an 
additional term, corresponding to an extra external field h 2 '. 



(ij) 



(22) 



The last two terms in the sum are introduced for further convenience. 

For the (n,m) HS, due to its self-similarity, after the summation over variables at intermediate sites we have the 
following equality: 



Z q (K,h,N) =exp 



ran 



Zll [K',h',— 
mn 



(23) 



where h = (hi, h 2 ). Expressions for the parameters of the new Hamiltonian K' and hf, and the function f q ^ may be 
found using the transfer matrix formalism. Introducing the transfer matrix for a bond: 



T mV2 = cxp[- J ff(ryi,77 2 )] , 
H {num) =K(1- 8 mm ) + y (2- 5 mfi - 5r, 2i0 ) 



Sr 11 ,oS rl2 ,0 - S 7hV2 - - (<5r fl ,0 + ^r )2 ,o) + 1 



(24) 



it is easy to calculate the transfer matrix for the (n, m) bundle: 



T' 





m 


(f n ) 


= exp 


\ / ?7lJ?2_ 





(25) 



with H' having the same structure as H, but with new parameters K' and h' . Then we come to the following equation 
for the free energy of the Potts model: 



f q (K,h) = ~lnZ q 



^-f q (K',h')-—f^(K,h), 



and for the percolation model: 



f(c,h) = —f(c',h')-—u(c,h), 
mn mn 



(26) 



where the variable c = exp(— K) is introduced instead of K, and u(c, h) — dfq ^ /dq\ q= \. 
Finally, after transition to infinitesimal transformations: 



, dX dX 

n = 1 + — , m=l + e— , 

A A 

dX 



dX 



dX 



c' = c + v c — , h' 12 = fti )2 +«i,2-r- , u = w— , 



(27) 



we arrive at the following equation: 



df df df 
oc ohx oh 2 



(28) 
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whe re v c , v\. 2 and w are found in Appendix |X| Eqs.(A4-A7). Setting h\ t % = in Eq.(28), and taking into account 



Eq.(A7), we have: 



X^r-(l + e)f = w =c+(l- c) ln(l - c) 



(29) 



where the independent variable A is related to c by Eq. (|l3|). 

Equation ( p8| ) should be supplemented by the boundary conditions. Directly from the definition ( |l9| ) we have: 
/ = — latc = and / = at c = f . Both Equations ( f28| , p9|) have essentially two different solutions depending on 
c > c t or c < ct, where ct is the threshold concentration of broken bonds determined from Eq. (14). 



IV. PERCOLATION EXPONENTS 



In the present section the critical exponents of connectivity will be obtained. Critical exponent v of the correlation 
length £: 



foe \c-ct\ u , 
is determined by linearization of Eq. (|l^) near c = c t : 

v~ x = e (1 + lnc) + 1 + ln(l - c). 
Eqs.([u|2|) can be solved analytically for eCl. In this case we have from Eq.(h 

After substituting Eq. ([52]) into Eq. (|3l|), the critical exponent v reads: 

1 

v w — . 

e 



(30) 
(31) 
(32) 

(33) 



The procedure for the solution of Eqs. ([l3 29) is described in Appendix^. The "partition function" /(c) is found to 
be: 



m = 



-c, c> c t ; 

-|(21n^) Hl 7(-i 7 l,21nt)- e i+i(lnt)'-, I ,-,,: 

■i-l,2Ina), c<c t ; 



H(21nt) i+1 r i 



(34) 



where T(a, x) is the incomplete T-function |40|, 7(0, x) — T(a) — T(a, x). It is curious that the "singular" part of the 
"free energy" to the left of the percolation threshold, c < c t : 



^ir (2e)' ■ - ( c t 

Is (c = 7T c t ■ * In — 

2 sm — V c 



7+1 



(35) 



is strongly oscillating function of e at e -C 1. The critical exponent for the "specific heat" a, f s cx \c — | 2 a appears 
to be large negative (as usual in the percolation model, see, e.g. ||). 



Differentiating Eq.(f28[) with respect to h% t 2 and setting hi t 2 = (see also Eqs.(A4-A7)) we get equations for the 
"order parameters" Pi (Eq.©) and P 2 = - df/dh 2 \ h=0 : 



X 



2-c+-(l-c) ln(l-c) 



cPi - P 2 



dP 

\^ = [c + \n{l-c))P x 



1 - c + -(2- c) ln(l - c) 



Pi 



(36) 



Boundary conditions: mi\ c=0 — 1, 777-2 1 c= o = 0' m i| c =i = m 2| c= i = may be obtained directly from the definition 
of order parameters. At c > c t we have the trivial solution of Eqs. (|36|), P\ — P 2 = 0. The lowest eigenvalue of 
the matrix in the right hand side of Eq.(|36|) at c = ct, which also is /3/V, may be easily found numerically for any 
dimensionality. The results are presented in Table | together with ones obtained by other methods. 
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For the nearly-ld system one can find the explicit form of order parameter dependence Pi (c) and Pi (A). Because 
the region of our interest is c < c t <C 1, Eqs. ( |36| ) can be rewritten for small c as: 



x dPl - c % c p 



A- 



dP 2 



-Pi + P 2 



One can see from Eq. (|38|) , that P 2 — 0(c 2 Pi), therefore one can neglect the second term in Eq. 
( |37| ) directly is solved by the substitution c = c t exp(— A c ) (see Eq. (Bl)) to read: 



(37) 
(38) 

As a result Eq. 



Pi = exp 



dx 

3e 1 ~ e 



-2x 



exp 



^Ei(-2A e ) 



where Ei is the integral exponent function. Taking into account the asymptotics of Ei (— x) = In a; — C 
x, and the relation between c and A, near the percolation threshold we obtain: 



(39) 



at small 



Pi 



c 2 1 

P = 7T = T" exp 
3e 3e 



(40) 



Thus the critical exponent (3 of the infinite cluster capacity appears to be very small and a strongly nonanalytic 
function of e. 

Using the scaling relations [|| , the complete set of critical exponents for the connectivity problems may be expressed 
through v and already obtained. 



V. CONDUCTIVITY EXPONENT 



Let us consider in more details the properties of the evolution functional B in Eq. (fill) ■ If we assume, e.g., that the 
conductivity is 0, with probability c, and equals some finite value (say, 1), with probability 1 — c, then its distribution 
function reads: c + (1 — c) e~ v . Choosing the distribution in Eq. ( jll| ) close to this form: 



P(v) = c+ (1 - c)e 11 +ip(rj) , 
where ^(77) is the small correction. Linearizing the evolution operator B with respect to tfi, we have: 
B ({P} , 77) = e (1 - c) Tie"" + e [c + (1 - c) e""] In [c + (1 - c) e""] + 



(41) 



d 2 ip 



dip 



(1 - c) In (1 - c) (1 - e-") + 77— \ + (1 - e) 77-^ + {l + e + eln [c + (1 - c) e""] } ^ . 

d?7^ d-q 

Another approximation for B is possible if we assume: 

P{rj) = c+{l-c) exp [-0(77)] , 



(42) 



(43) 



with </> (0) = 0, (f> (77) — * ±cxd as 77 — > ±00 rapidly enough (faster than ±a/|77|, as we shall see later). 

The important point is also to assume analyticity of P (77) and of Q(s) within some stripe along the real axis. Using 
the relations: 



r( 2 ) 



1 ' ) (z) = as z — > , 

7TZ 



(12) -. 

where i?^ ' are Hankel's function of first and second kind, respectively, we may replace the integrals with the 



Ji-functions in Eqs. (H,|ll]) with the ones containing Hj 1 , along the contour C shown on Fig. [| Thus we obtain: 



Q(s) = (1 - c) 



/>00 T 

Vi/ -^:Ji(2V^7) exp (-0(77)) 
Jo 

= -{l-c)^J c (2^j) exp (-^ (77)) , 



(44) 
(45) 







FIG. 3. Integration contour in Eq. 

where the contribution of the pole of (2y/sTp) /y/rj in Eq. ( |45| ) just reproduces the first term in the square 
brackets of Eq. (|44"|). Assuming \s\ to be large enough, one may replace by its asymptotic expression: 

H[ 1} (2^) * k- 1/2 ( S77 )- 1/4 exp (-^ + 2i 

and treat the integral (Eq) by the saddle point method. As a result, we get: 



Q(*) = (l-c) 



4>'(vc) 



-,1/2 



exp [-0(??c) + 2?7 C 0'(77 C )] 



>'(r ?c )+27 ?c ^"(? 7c ) 
where the saddle point should be determined from the equation: 

i\—-4>'(Vc)=0, or ?7c (0'(f? c )) 2 = -s . 

V 

If we define: 

X ( s ) = fac) - ZiyfsTh = (??c) - 2r? G 0'(rj c ) , 
the following relations can be easily established: 

1 



(46) 



(47) 



(48) 



1 

x'(s) + 2sx"(s) = 



(j)'(r] c ) + 2i lc (j3"(7 lc 



Obviously, the transformation 4>(j]) <-> x(s) is symmetric, i.e., its reverse has the same functional form. With x(s) we 
may rewrite Eq.ffiq) as: 



Q(a) = (1 - c) 



l + 2.s 



-,1/2 



exp(-x(s)). 



(49) 



Proceeding along the same line, one also can derive the following equality: 

y/tj / -=J X (2y^) Q( s ) lnQ(s) = (1 - c) In (1 - c) + 
Jo V s 



(1-c) 
(1-c) 



x(Sc )-ln ( l- C )-iln(l + 2 Sc ^) 



exp(-x(s c ) + 2s c x'(s c )) 



(50) 



0(r ? )-2#'( 7? ) + -ln I + 25- 



'(»?) 



'(V) 



exp (-^(77)) 
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Replacing Mrj) = — ln[P(r?) — c] back in Eq. (|5C|), and substituting the result into Eq. (|ll|), we get the evolution 
equation (y_lj) in the differential form: 



dP 

A— = B x ({P} ,rj) = - (1 + e) r?P' + ePlnP - (1 - c) In (1 - c) + (P - c) In (P - c) 



(P - c) In 



P" P' 
1 - 2 V — - 2 V 



P-c 



(51) 



The derivation of Eq. (|5l]) suggests the replacement of the evolution operator B with its approximate form B\ 
at least for large enough r\. But, if the function P(r\) is analytic, then, taking into account P(0) = 1, it may be 
represented in the form ( |4l| ) in some neighborhood of the point r\ = 0. If we plug Eq. ( f4l| ) into Eq. ( |51), a nd linearize 
the resulting expression with respect to ip, we arrive exactly at the same evolution operator as in Eq. (42), which was 
obtained by the linearization of the exact evolution operator B. This observation prompts us to enhance the region 
of validity of the approximate evolution equation ( pl[ ) for everywhere in the complex plane rj. 

At the percolation threshold, c — c t , the solution of the RG evolution equation can be taken in the form: 



P(r),\) = P{riX~ a ). 



(52) 



Here the critical index a is related to the critical exponents of the conductivity fi and of the correlation length v by 
the relation 

a = [ijv. (53) 
Indeed, according to Eqs. (0,p3), the average conductivity of the A-sized cube in the critical regime, c — Ct, is 



< <r(A) >= 



rfPfa A) 
dr) 



A 

|,,=o — c \ — 



(54) 



The same conductivity (Q) is realized in the infinite disordered fractal with the correlation length £ equal to A. Near 
the percolation threshold the correlation length £ is given by Eq. (^o|) and, therefore, according to Eq. (|54|) the fractal 
conductivity obeys the scaling law 



fi = av. 



With the scaling anzats (|5^), Eq. ( pl| ) becomes an ordinary differential one of the second order. It appears to be more 
convenient to use the function <f)(x) = — In [(P(x) — c t ) /(l — c t )] instead of P(x). Denoting <fr — — In [c t /(l — c t )], 
we have: 



In 



1 + 2x - 



where we introduce: 



(1 + e - a) x<f)' - 4> + e [g (0 - <fo) - g (-</> )] 



g(<t>) = (e* + 1) In (l + e*) 



(55) 



(56) 



An equation for <fio which follows from Eq. ( |l4| ) was used to derive Eq. (55). The latter may be reduced to the first 
order linear equation by introducing: 



z (0) = cxp [-2 (1 + e - a) x<f>'] , 

and treating (f> as independent variable: 
dz 

— + (1 + e - a) z = - (1 + e - a) exp {-20 + 2e [g {cf> - <p ) - g (-&)]} 



(57) 



(58) 



where a is related to the critical index of conductivity by Eq. (p3|). 

One should require the solution of Eq. (58) z((f>) — > as <fi — > +oo faster than exp [— (1 + e — a) 0] to ensure the 
applicability of the saddle-point approximation. This selects a solution in the form 



/>oo 

z (cj>) = (1 + e - a) e-^- a ^ / dy exp {- (1 + e - a) v + 2e [g (y - O ) - 5 (-<^o)]} • 



(59) 
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The normalization condition 0(0) = implies z(0) — 1, from which the equation for a follows: 



(1 + e-a) dyexp{-(l + e-a)y + 2e[g(y-<f> )-g(-<f> )}} = l. (60) 
Jo 

Comparing the values of a, obtained from Eq. (|60|) and by the numerical investigation of the original evolution 
equation (pi]), one can see that both methods give the same results at any dimensionality. This, together with the 
considerations presented above, prompts us to consider the saddle point solution as an exact one. Of course, the 
RGMK method itself is an approximate one. 

Comparison of the numeric results for the critical exponent of conductivity /i is presented in Table |. In three 
dimensions we have from Eq. (60): a w 1.891. On the other hand, the best numerical results ^ij give a = 2.25±0.04. 
So, the RGMK method provides a reasonable solution even for 3d systems. 

For the case e«l, from Eq. (|6C| ) it follows 

1 + e ( l + e \ a fm\ 

a = cxp , fi = —. (61) 

e \ e J e 

For e ^> 1 from Eq. (|6(]) one may obtain: 

a = e- ( -e^. (62) 



VI. DISTRIBUTION FUNCTION AT THE THRESHOLD 



The function 



may be determined as a reverse of the equation: 



Cx 



>exp 



<K 



*(0 



(63) 



The integration constant C in Eq. flS3) corresponds to arbitrary choice of the unit of conductivity, or, alternatively, 
of the length scale at the percolation threshold. 

Thus the distribution function (DF) for conductivities which in the initial representation is defined as, n(cr, A) = 
(S (a\ — a)) at the percolation threshold takes the universal scaling form: 



n(a,\)=ct + (l-c t )n(y), 



where y is the conductivity in units of the average conductivity 



< cr(A) > 



a 

00 



and the scaling function Tl(y) may be expressed as the integral: 



n(y) 



dx 



1 



2m 



V Jo 



0+2OO 



— exp [xy - 4>(x)] = - — exp [-</> + yx(<f>)] 



2ni 



(64) 



(65) 



(66) 



The last equality was obtained through integration by parts. 

However, it should be noted that additional unphysical contribution arises when evaluating the integral in Eq. (|6f 
Namely, the function x((j>) is singular at cf) — </>„ = <fio + iir(2n + 1), where n is an integer number. These singularities 
result from the procedure of analytic continuation within the RGMK approach. It can be illustrated as follows: 
Let us assume the initial distribution of conductivities to be: Po(fl) = c + (1 — c)e~ n . After putting m identically 
distributed conductors in parallel, the Laplace transform of DF for their sum, Pi(r]) = [c+ (1 — c)e~ v ] m , has m-th 
order zeroes at fj n = r)o + iix(2n + 1), where 770 = hi[(l — c)/c]. These zeroes transform into singularities after analytic 
continuation to non-integer to. Thus the procedure of the transition from integer rescaling transformation (which is 
exact for a hierarchical structure) to the infinitesimal one is the source of the above singularities in Eq. (pw. Since 
these singularities are artificial ones they should be merely discarded in the integral (|6^). 
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At large conductivities y 1, shifting integration contour in Eq. ( pq ) to the region > fa, one has the following 
asymptotic expression for the DF: 



V2 \V2 



exp 



where: 



(67) 



[2 7 r(l + e)(l + e-a)] 1 / 2 fl + e + a\^+ 
(1 - c t )a 



A 2 = 2{l + e-a) f dCM~0 

•/ — oo 



1 + e- a 

j C 
dC lnz(C)' 



, 2/2 



1 + e /l + e-a\ 1+£ a 2 

e ^ 



1 + e- a 



(68) 



Shifting the integration contour in Eq.(|66j) to the region 3t(f> < fa, we arrive at the following expression for H(y) in 
the region of small y: 



D x ( y 



U(y « 1) = -- l± 

vi \yi 



7 + 1 



exp 



(69) 



with: 



e — a 



Di= IM1 + e _ a) |V % a)| 

e — a 1 + e — a 11 + 6- a 



/•OO 

Ax =2(l + e- a) / dClnC 
Jo 



j c 

d( lnz(C) 



(70) 



More detailed results are available in the limit e <C 1. At 3t<p < fa ~ 1/e to first order of Ct and of a we get the 
following expression for x{(j>): 



\rxx{4>) = In (f> + ct 



e*-l- 



j(i+«)C _ i _ (1 + e )C 



(71) 



Evaluating Taylor's series of (a;) at x = 0, the central momenta of conductivity are found to be of the order of a: 



(*-(*))' 



a + c t . 



= _ t:(1 + e ) a + c t ■ 



(72) 



On the other hand, using in Eq. (|66j) the asymptotics of x{fa) at $t(j> < fa and \<f>\ 3> 1 we have for large enough y 
(see also Appendix : 



where 



yU(y) = e^^^-WSiW), 
a 

c t , and the new fluctuating variable was introduced: 



(73) 



1 + e a 



7 is Euler's constant, and S(W) is given by: 

S(W) 



zoo+A 



ioo+A 



'i 

2ttz 



(74) 



(75) 



An asymptotic expression for S(W) may be easily obtained by the saddle-point method (Appendix |^) 



S(W) 



cxp(-e~ 1 VK) 
V2ireW ' 



as W » 1 : 



In ln(e/Vy) 
ln(e/W) 



, as W < 1. 



(76) 
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Fig. @ shows WS(W) as a function of ]nW. 



w s 

. 2 r 




In W 



FIG. 4. Universal distribution of conductivity for the quasi-lD fractal at the percolation threshold (Eq. (^)). There is a 
sharp decay in the region of large conductivity and a long tail for large resistances (Eq. (|7rj)). 



The asymptote of n(y) at small y is given by Eqs. ( p9| , |70| ). At e <C 1 it may be reduced to: 



yfl(y) 



exp I - - e x y l > ' 
2ire \ ^ 



(77) 



The two expressions ([73|,[77|) should be supplemented by one for the intermediate region, where the function x(<j)) in 
the integral (pq) can be expanded in a and Ct- Here we have in the first order of a, ct- 



x(cf>) = + r? (0) , r)(4>) =ct(e*-l-0) + 



(i + e ) 2 7o c 2 



o(l+<0C 



l-(l + e)C 



yn(y) =S(y-l) 



2tti 



Apart from the (5-function term, this yields in the region < y < 1 : 

a 1 



n(y) 



1 + e A 2 ' 



(78) 



where A = 1 — y. 

To establish regions of validity for different expressions of DF, let us consider the region Ac 1. Here Eq.(f77|) turns 
into: 



1 



1 A 



exp — e 

2ire \2 2e 



A/e-l 



n(y) 



and Eq. (|73|) may be written, taking into account Eq.(|76[) at W <C 1, as: 



n(y) 



/27T a 
e 1 + e 



In (In ±±£ + A) 
A + ^ln% 



(79) 



(80) 
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Comparing Eq.@ with Eqs.@^) one can conclude that Eq.@ is valid if A < Ai ~ a/eln(l/e), and Eq.@ 
holds for A > A2 ~ e In (e /e) . Fluctuations of conductivity appears to be distributed within narrow region of relative 
width Ai, which ensures that not very high order central momenta of the conductivity to be small (see Eq. (|72|) ) . 
However, if expressed in terms of the universally fluctuating variable W , the distribution becomes smeared over a 
wide region with the lower cut-off W\ ~a p ,j)»il+l/ln (V 6 )- 

The distribution function S (W) arises naturally in a Id chain of random resistors, if, to require a scaling form for 
the distribution function of A-length chain: T (p, A) = Y(pX~ a ), or Q (s, A) = Q (s\ a ) in the Laplace representation. 
Then from Q (s, nX) = Q n (s/n, A) it immediately follows: Q (x) = exp [— Cx 1 /( 1+a '] . Evaluating its inverse Laplace's 
transform T (r), and assuming a <C 1, which is true in the Id case, we have after the proper rescaling of the integration 
variable: 

r T(r) = -WS(W), W^ar- 1 ^, (81) 
a 

which is essentially the same formula as Eqs. @|74]). 



VII. SCALING AND AC CONDUCTIVITY 



Exact results for AC-conductivity in disordered systems are available for a very limited class of models, mostly for 
Id ones. The common method to study disordered hopping systems is the effective medium approximation (EMA), 
which gives qualitatively correct results for three-, two-, and even for one-dimensional systems. However, it fails for 
a nearly Id system. In a percolation model, for an example, EMA gives threshold concentration value Ct oc e and 
completely wrong values of critical exponents. However, knowing the results for DC conductivity and topological 
properties of the percolation network, the qualitative behavior of low-frequency conductivity may be restored within 
the scaling hypothesis Hfl2) ■ 

Namely, it should be assumed that the only length scale near the threshold is the correlation length £ cx |t| _i/ , 
t = (ct — c) I Ct ■ The second assumption is about the anomalous diffusion of a tracer placed onto the infinite cluster 
at the percolation threshold: 

(r^^oc^, C<1- (82) 

r(t) is the distance from tracer's position at t = 0, (. . .) means the average over the initial positions within the 
infinite cluster only. Above the threshold, when c < c t , we have normal diffusion at sufficiently large times, when 
(r 2 (t))>e: 

(r 2 {t))=D OQ t, (83) 
with the diffusion constant D^, connected with DC conductivity as: 

aDC = e -w Da °- (84) 

e and n e are charge and concentration of electrons, respectively. Near threshold we have cx r M , /1 is the critical 
exponent of the DC conductivity. Note, however, that in Eq. J83] ) the average is over the whole network, including 
finite clusters. 

The relationship between exponents /i and £ may be established in the following way: At some large enough time t 
let us consider the A-sized box, A = <V 2 (£))^ 2 oc t^l" 1 inside the system at the percolation threshold. Its conductivity is 
tr(A) oc \~^l v . At the same time, it may be expressed as tr(A) oc P OD (X)X 2 /t, where Poo (A) = Pi is the infinite cluster 
capacity for the system in which the correlation length equals A. Since (A) oc A-^, it follows o-(A) oc X^+^K . 
Thus we can conclude: 

2v + p — p 

If we take into account the average over finite clusters, we obtain at the percolation threshold: (r 2 {t)^ oc Poo(A)A 2 oc 
^2-/3/1/ (x ^ i s the anomalous diffusion exponent including the contribution of finite clusters: 

( .( 1 -') ( -»dL- (86 , 

V 2v J 2v + p — (3 
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Below the percolation threshold, at c > Ct, we have (r 2 (t)') ~ ^ 2 a t 2v at t — » oo. 
All the above may be summarized as: 

(r 2 (t))=^G(Tt u ^), (87) 

where: 

« = 1 - C = 9 / „ , (88) 
2v + [i — p 

with the scaling function G(0) being some constant, G(x) cx x^ as x — > +oo, and G(x) oc | a? | s , 

s = 2i> - /3 , (89) 

as i -> — oo. 

The conductivity may be expressed through (r 2 (5 1 )), the Laplace transform of (r 2 (i)), as: 

*(u,,T) = ^S*(r 2 (S)) , (90) 

where 5 = —iu. Using the Tauberian theorem for the Laplace transformation of power laws, from Eqs. ( |87| , |9"o| ) one 
can obtain: 

*{v,t) = j£-S»g(tS-«''>) , (91) 

with scaling function G(x) having the same asymptotic properties as G{x). 
Thus we have at c < ct and u> <C r s+fl : 

at c> c t and u> < |r| s+M : 

(J CX — iU!T~ S , 

and at cj 3> |t| s+ ' 1 (in particular, at any a; if c = q): 

ct cx (-iu>) u . 

The summary of the frequency dependence of AC conductivity is given by Fig. [5]. 
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FIG. 5. Diagram of frequency-dependent conductivity in the quasi-lD fractal near the percolation threshold, |r| <C 1. Here 
s = 2v — (3 and u — fi/(2v + fi — j3) and v, (5 and fi axe the critical indexes of the correlation length, capacity of infinite cluster 
and conductivity; ui c ~ |rj s+M is the boundary frequency separated the critical region from the conducting and dielectric phases. 
For the quasi- ID fractal with the transverse dimensionality £<l,n<l,|ti<l and s> 1. 



In a nearly Id case the static conductivity exponent [i is given by Eq. ([H]), the frequency dependence exponent 
and the exponent s of dielectric constant divergence may be written with Eqs. ( ]33| , [40| ) as: 



u 



2 

s » - 
e 



(92) 
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uw^exp^-l-^ . (93) 

Thus the exponent s is very large but the exponent u is very small. Therefore, in a dielectric phase the AC conductivity 
as a function of frequency demonstrates a step-type behavior. In the conducting state the frequency dependence of 
conductivity remains very weak. 



VIII. VARIABLE RANGE HOPPING 



In a conductor with localized carriers the charge transport is provided by the variable range hopping (VRH). The 
model may be formulated as follows. The phonon assisted hopping rate my from one localized state j to the other i 
per unit time, including Fermi occupation probabilities pi , is approximated by the formula |17[ : 



jPj (! ~ Pi) = w exp (-2/y) , 



Hi — 



£i - £f\ 



AkT 



(94) 



where e, and are energies and position vectors of localized states respectively, and a is their radius. We assume here 
the hopping motion to be along the chains. Assuming that the localized states near the Fermi level are distributed 
uniformly in space and energy, the distribution of random variables /y is: 



F(f) = Probability^ > /) = exp 



where: 



fo 



1/2 



kT 



1 



AN F a 



(95) 



(96) 



Np is the density of states at the Fermi level. 

To study charge motion in a system with continuous distribution of hopping rates is a much more complicated 
problem than the one for a percolating system, where hopping rates are either 0, or some given finite value w$. 
However, knowing the results for the conductivity of percolating system, qualitative conclusions can be obtained 
for the system with continuously distributed hopping rates. Namely, let us introduce some probe hopping rate w c , 
replacing all hopping rates my < w c with 0, and all my > w c with w c . Obviously the conductivity becomes lower 
than the initial one, but if to choose w c from the requirement to get a maximal conductivity, one can hope to obtain 
a good estimate for the original conductivity. 

It is convenient to represent the probe value as w c — exp(— 2/ c ). The corresponding broken bonds concentration c 

is given by formula (|95|), i.e., c = exp — {fc/ fo) 



The value of f c for the threshold concentration is f t = /o/ye or 



ft = 



1 / T\ 



1/2 



4kT Q 



2 V T J e eN F a 

Assuming the probe value f c = f t + 8 to be close to the threshold one, r = (ct — c) /ct <C 1, we have the relation 

s=y t r. 

The scaling formula ( |9l|) for the conductivity of the percolation system now reads 



(97) 



cr(w,r) 



e 2 1 iji 

—a\\\rfw c g 



kT 



(99) 



where w c = w t exp (— e/tr), wt = wo exp (— 2f t ), aii = l/{NpkT) is the hopping length, and the electron density 
n e = NpkT. Scaling function g(x) has the following properties: 
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A I 



as \x\ < 1 , 



D + B + x- s -f* 



. . . as x >■ 1 , 

as a; < 0, |a;| ^> 1, 



with coefficients A, B± and Z? of the order of unity. 

First let us consider DC conductivity. Obviously one should choose r > 0, and the expression to maximize the 
conductivity as a function of r is t m exp (— e/tr). The optimal value of probe parameter r is very small: t^c = eft- 
Taking into account Eq. (|9^), we have: 



0~DC 



7V F (fcT) 



rWt, 



lUi = Wo ex P 



1/2' 



(100) 



Thus the DC conductivity obeys a quasi-ld Mott's law, but the characteristic temperature T\ given by Eq. ( |97| ) is 
much greater than To for VRH in a strictly Id chain from Eq. (p6|). 

The so-called "hydrodynamic region" of very low frequencies, where the conductivity's frequency dependence is 
determined by expansion: 



<y = <jdc i 

UJh 



UJ <UJ h , 



(101) 



appears to be very narrow for the nearly Id system. Its width u>h may be estimated from the condition on the 
argument of the scaling function g in Eq. (|99|) to be of the order of unity at r = tdc and ui = u>h- 



2/6 

T DC W t 



exp 



,-2/e 
ft w t- 



(102) 



Within this region \uj\ < u>h we get the effective value of r c to be dependent on frequency as r c = tdc (1 — S/u>h + ...). 
We suppose at further derivations S — —iu> to be real and positive, having in mind analytic continuation afterwards. 

At S ~ u)h, t c changes its sign, and now the conductivity is determined by the charge motion inside finite-size 
clusters. At |r c | <C 1, the size of effective clusters is large, i.e., the clusters contain many Id chains. This frequency 
region is called multiple hopping one. Note, that in contrast to two- and three-dimensional systems for which the 
multiple hopping regime transforms at higher frequencies into the regime of pair hops, here the multiple hopping 
frequency region borders that of the one-dimensional hopping. 

From properties of the scaling function g in Eq. (|9{]) one can conclude, that the conductivity is maximal if one 
chooses the probe value r c such that the argument of g is of the order of unity. Taking into account explicit expressions 
for critical exponents in a nearly Id case we have: 



r exp 



e 2 ftr' 



(103) 



where r' = — r c and the conductivity may be estimated to be: 



Np(kT)'- 



C T '-2/e 



N F {kT) 



r wtexp(e/tr') . 



(104) 



One can see from Eq. ( |103[ ), that the character of frequency dependence is determined by the parameter e 2 / t = 
2(T 2 /T) 1 / 2 , where: 



T 2 = -e 4 T! 



(105) 



If the temperature is relatively high, T S> T%, from Eq. ( |103| ) it follows r' ~ {S/wtY^ 2 , and the conductivity reads: 

,,2 



a 



N F (kT) 



2 w t exp 



5W 2- 

£ft[ m 



(106) 



Thus the temperature dependence of conductivity at a given frequency is described by the quasi-ld Mott's law 

1 /2 A 1 

exp — (Ti/T) 1 . The frequency region for application of Eq. ( jlOq ) is determined by the requirement r < 1, which 
may be written as: 
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5 h 2 »> 



At lower temperatures, T <C T 2 , Eq. (106) is valid as long as (l/2)e 2 / t r' <C 1, or: 



6 1 w * ^ 1 6 1 1 T2 

_ ln _» ln _ = _ ln _. 



At higher frequencies the solution of Eq. (102) is given by the equation 

, 1 , S 



r — In — , 

eft Wi 



2 /t 

where uj\ — (2/e 2 ft) = (T ' jT^ 11 Wt, and the conductivity now is 

\ 2/e 



2A 




N F {kT) 2 \>£ 
This formula remains to be valid until r' <C 1, i.e. if: 

- In — < eft • 

2 w t 



(107) 



(108) 



(109) 



(110) 



If the frequency is higher than ones determined by Eqs. (107) or ( |11C| ), the conductivity behavior becomes a Id 
one (see, e.g., Ref. fllif ). 



IX. CONCLUSIONS. 



As one naturally expects, the results for the percolation problem in nearly one dimension approach Id ones as 
e — > 0. In particular, threshold concentration c* (Eq. fl32|)) and critical exponents (3 and /i (Eqs. ([l0 61)) tend to 



zero. As a result the capacity and conductivity have a jump-like behavior as a function of concentration of broken 
bonds near a critical value. The reason is that in the limit of one dimension the infinite cluster arises at c = and 
occupies immediately the whole system. The critical length exponent v « e _1 is, however, large, contrary to the Id 
case, when v — 1, but this nearly Id behavior of the correlation length can be observed in a very narrow range of 
concentrations, |c — Ct\ <C 1. Outside this region, when C(«c«l, critical length scales in a Id manner, £ = c _1 . 

The other surprising feature is the strongly nonanalytic behavior of both threshold concentration and of critical 
exponents, which points to a regular e-expansion near lower critical dimensionality d = 1 being rather impossible. 
Although the RGMK method becomes exact only in the limit e <C 1 the comparison with numerical results (see Table 
1) points out that the critical indexes obtained by this method proves realistic even for e = 1, 2. 

All these features together teach us that the infinite cluster arises almost like a jump. The infinite cluster at the 
percolation threshold itself is a fractal with a number of dimensionalities B all of them less than the dimensionality 
of the original lattice itself. For example, fractal dimensionality: 



D f = D 



t 
v 



1 



1 



exp 



(111) 



characterizing the mass distribution within infinite cluster, is very close to the fractal dimension of the system itself 
D, which means that infinite cluster at the threshold is "almost dense". 
Fracton, or spectral dimension d: 



D 



2v 



2v + fi - (3 



1 + e 



1 

2~c 



exp 



■1 



(112) 



was introduced to describe the behavior of random walk on the infinite cluster (it may also be used to describe, e.g. 
density of localized vibrational states, or fractons, etc.). Its closeness to D = 1 + e means that the diffusion on the 
infinite cluster at the threshold is almost normal. Respectively, the conductivity frequency dependence exponent it, 
Eq. (|93|), is small. On the other hand, dielectric constant in the insulating phase, e' cx \t\~ s diverges very strongly, 
nearly as £ 2 oc |r|^ 2ly , but, similar to the correlation length £, this divergency takes place in a narrow interval of 
concentrations of the order of Cj. 
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The RGMK enables us to study not only the average characteristics of the system but also their fluctuations. These 
fluctuations become essential near the critical point when the correlation length £ becomes larger or comparable 
with the system size A. In this case because of non-self-averaging a sample demonstrates individual characteristics 
corresponding to its specific disorder. 

In the present work we have found the distribution of possible conductivities of samples in the critical regime. The 
average conductivity < er(A) > decays with the sample size A according to scaling law All the fluctuations are 

found to obey the same scaling law ([72]). Thus the distribution of conductivity is the universal function, fl(y), in units 
of the average conductivity, y = a j < <r(A) >. In other words the above fractal dimensionalities of the percolating 
cluster do not vary with the fractal size, A, as happens in multi-fractal systems f43] |. This robustness comes from the 
additive laws (|J) for classical charge transport. 

The function tl(y) represents the distribution of possible experimental deviations from the scaling law (p>4|). It is 
shown, Eq. ([72]), that the central body of the distribution Tl(y) is concentrated in the narrow interval around the 
average value, i.e., y = 1, but it does not take the gaussian form. The distant tails of the distribution in the region 
of large conductivity y » 1 and for large resistivities l/y> 1 decay like a stretched exponent, Eqs. (|67|,|69|). In the 



limit e«1 the shape of the distribution function (|73[) is consistent with the Id scaling of the percolating cluster. 

Returning to the variable range hopping in the chain fractal as a consequence of (i) Id character of variable range 
hopping along the chain s, an d (ii) finite value of broken bonds threshold concentration c t , the DC conductivity obeys 



a quasi-ld Mott's law (100). But the characteristic temperature T\ of this dependence is higher than the formal 
value of the characteristic temperature for Id chain (remembering that Mott's law is not valid for Id systems), T Q , 
by a factor 1/e. This increase can be understood through comparison with the quasi-ld model of weakly coupled 



metallic chains 14 1. The variable range hopping conductivity of this model obeys the same law with the characteristic 
temperature T* = 2 (d-i) ■ wnere 2(d— 1) is the number of neighboring chains and d — 1 is the transverse dimensionality 
of the quasi-ld system (d = 3). Taking d — 1 = e we formally reproduce T\ for the nearly-ld fractal. Experimental 
temperature dependence of conductivity in the poorly conducting polymers very often follows a quasi-ld Mott's law 
with substantially increased characteristic temperature [ ^fj|Jll|j44]| , as expected for eCl. 

In 2d and 3d isotropic systems with VRH mechanism of charge transport the temperature dependencies of the 
AC conductivity and of the dielectric constant are rather weak, but in nearly Id systems there exists the region of 
f requ encies and temperatures, where these dependencies are nearly the same as the quasi-ld Mott's type (see Eq. 



( |106| )), which continuously transforms into Id dependence within a wide enough transient region, Eq. (109). Such a 
type of strong temperature dependence for both DC and AC conductivity and also dielectric constant is experimentally 
observed in conducting polymers with localized carriers |^,[^] . 

The physical picture behind this dependence is the following. The low-dimensional random system can be separated 
in weakly coupled clusters within which carriers are confined. With increasing temperature, the size of clusters 
exponentially increases, as more space accessible for carriers due to thermal activation. As a result, the dielectric 
constant and the conductivity exhibit strong temperature dependencies. In contrast to the low-dimensional case, the 
clusters in two- and three-dimensional systems prove to be more effectively coupled. Therefore the large polarization 
of clusters does not happen because of transition of carriers between clusters. Thus our results support strongly the 
idea that even poorly conducting polymers represent low dimensional systems. 
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APPENDIX A: EVALUATION OF COEFFICIENTS IN THE RG EQUATION FOR "FREE ENERGY". 

Comparing Eqs. (??) and (p5h, we have: 



4°)(c,/ l ) = ln^ 0! u(c,h)=^n 



rp/2 

c! = e- K ' = -4^ 



q=Q * 00*11' 

2 



rpl rpl / rpl * 

J oo i oo \ 1 wJ 

where i, /' = 1, . . . , q — 1, and the results are independent of the choice of I ^ V . From Eq. (??) it follows that: 
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(f n+1 ) = T 00 (f n ) +(q- 1)T Z0 (f n ) = (f n ) +(q- l)ce- (h ^ h2 ^ 2 (f n ) , 
V /oo V /oo V J 10 V /oo ~ V /;o 

and, analogously, 

(t" +1 ) = ce-<- hl - ha V 3 (r") + e~ hl [l + (?~ 2)ce /l2 ] (r n ) ;f) , 



(T ™) + e^ 1 f f n ) + (g - 2)ce~ 



, — hi+h 2 



(f n+1 )^ = ce-(^-^)/ 2 (r") +ce-' ll+/l2 (f") H +e-' 11 [l + ( g -3)ce' 12 ] (t") . 



Here the symmetry property T m „ 2 = T„ 2711 was used. Introducing: 



u(n) = |- (f 

o<7 V /oo 



*i(n) = T 



-3=1 



9=1 



, i 3 (n)=(T" 



9=1 



and taking into account that at q — 1 we have T 00 = ( T n ) = Tn = 1 , the following equations can be established 

V / oo 

to read: 

h(n + 1) = ce -C»i-^)/2 + e -h, ^ _ ce h 2 j ti ( n j j 

fe(n + 1) = ce-^ 1 -^^ 2 *!^) + e- h n 2 {n) - ce- hl+h H 3 (n) , 

t 3 (n + 1) = ce-( /ll -' l2 )/ 2 ti(n) + ce-' ll+ ' l2 i 2 M + e~' 11 (l - 2ce' 12 ) t 3 (n) , (A2) 
u(n + 1) = u{n) + ce-( ftl_ ' l2 )/ 2 ti(n) , 

with initial conditions: ti(0) = ^(0) = u(0) = 0, ^(0) = 1, which can be easily solved. After the transition to 
infinitesimal transformation, n = 1 + d\/X, rn = 1 + edX/X, we have, taking into account Eq. (p5|): 



T'oi = t?(n) = ce -( ft i-M/ 2 J i + 



r/A 



e In c 



hi] 



T' u =^{n) = e- h ^l+^ 



c 2 e 2h! + h 2 



(e hl + ce h * - 1) 



-ehi + ce h2 ( 1 - 
2 J t ln (! 



e^ 1 + ce h2 — 1 



(A3) 



T' w =t^{n)=ce 



1 + -Tr- 



ee: 



e (lnc - /ii + fa) + 1 - ■ , 



„2/ii 



| lnT °° 



r/A 



c 2 e h2 



[m(i 



9=1 



A (e^i + ce' 12 - If 



1-1 



e ftl + ce h2 



fa 



From Eqs. (A3), (27) and (Al) it immediately follows that: 



v r = ec In c — c I 1 — 



e /ii + 



26^ 



hi _|_ ce /i 2 _ 1 



ce' 12 — 1 
• 1 



,2/ii 



[ln(l 



ce 



h 2 \ 



fa 



(A4) 



22 



Vi = ehi - ce hl 1 



,,/u 



,hi _|_ ce h 2 _ I 



c 2g2/ii+/i2 



[in (1 - ce h *) - hi] , 



(A5) 



v 2 = eh 2 + (2 - ce h2 ) ( 1 



ce 



hi 



c 2 e 2h 1 +h 2 



e hi _|_ ce h 2 _ i 

2e hl 



{ehi + ce h 2 _ if e fci + ce^ 



1 



"la _ 1 



[In (1 - ce 712 ) - fa] , 



(A6) 



(e^ + ce h2 - I)' 



1-1 



e™ 1 + ce'' 



__ [ m (l-ce>-)-fa] 



(A7) 



APPENDIX B: EVALUATION OF TOPOLOGICAL QUANTITIES AT e«1. 
At c <C 1 Eq. ( |l3| ) may be rewritten as: 



which can be easily solved to yield: 



dc 



c t exp 



ec In • 



(Bl) 



where Ao is the arbitrary positive integration constant. On the other hand, one can neglect first term in Eq.(|13[) when 
c»Cf, obtaining the solution: 



c = 1 — cxp(— A/Ai) , 



(B2) 



with the other integration constant Aj. 

One can match plus sign solution (|Bl| ) with Eq. (B2) in the region q < c < 1. Setting Ao = e 1 / e Ao in Eq. (p3l|), 
one can see, taking into account Eq. ( |32[) , that to fulfill f ^s> c ct, one should require < 1 — (A/Ao) c <C 1; therefore 
one can set (A/Ao) e ~ 1 + eln(A/Ao), and formula (Bl) turns into c w X/Xq. On the other hand, we have from Eq. 
(B2), that within the same region c w A/ Ai, that is, Ao ~ Ai, or Ao ~ e 1 / e Ai. 

When c«l, taking into account (Bl), Eq. ( [29|) reads (we set the integration constant Aq = 1): 



A 



dX 



(l + e)/ = yexp(±2A e ) 



(B3) 



At c < ct (minus sign in the right hand side of Eq. (B3)), one should set / = at A — > oo, and the corresponding 
solution is: 



/(A) = ^(2A e ) 



1 c\ (2A £ ) 

2 1 + esin^ T (■ 



7 + 1 



2 00 



(-2A e ) r 



n![l - e(n- 1)]' 



(B4) 



where r(a,ir) is the incomplete T-function |40[| , tending to at x — > oo. Substituting in the first term in the right 
hand side of formula (B4) A € = ln(c t /c), an d re placing T(l/e) with its asymptotic expression, one arrives at Eq. (|3F 
At 1 ^> c > ct, the general solution of Eq. (B3) with a plus sign is: 
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/(A) 



(-2A e )- +1 7 ( ---1,-2A C ] - AX 



2c 



(B5) 



where 7(0,2;) = T(a) — T(a,x) (note that x~ a j(a, x) is analytic function of x at x = 0), and A is an integration 
constant to be determin ed through combining with the expression for the "free energy" at c > c ( . The latter may be 
done, substituting Eq. ( JB2| ) with Ai = e 1 / 6 into right hand side of Eq.(p9|). 

Having in mind that solutions (B+2^) are to be matched with each other in lowest order on e, and setting e = in 
the left-hand side of Eq.(p9|) too, we arrive at: 



/ = exp ( -e 1/e A ) - 1 



(B6) 



On the other hand, from Eq.(B5) we have within the same region, replacing the function 7 with its asymptotics at a 
large value of its second argument Gof : 



4eA' 



exp (2A e ) - AX 



l+e 



4e In • 



Ae- 1 ^ 



(B7) 



Comparing expressions (B£) and B7), one can conclude that A w e 1 / 6 " 1 " 1 . Substituting this into (B5), replacing there 
A e with ln(c/ct), and combining with Eqs. (B4) and (B6), one arrives at Eq. (B4). 



APPENDIX C: INVERSE LAPLACE TRANSFORMATION OF THE DISTRIBUTION FUNCTION. 

After some integrations by parts, Eq.(|7l]) may be rewritten as: 

lnx (cb) = wr+ ( 1 — ] hi (-</>) H — [ 1 - 7 - ln(l + e ' 



l + e J v ' l + e V + 

(1 + ef .Loo C 2 

»t7r+(l-— - + Gi+c t - r + - 3— r— » C1 

V l + e/ l + e cb (l + e) 4> 

4>\ = a I (1 + e) 2 , Gi = 1 —7 — In (1 + e) — (1 + e) Ct/a. From the letter expression in Eq.( |Cl|) , which is the asymptotics 
as \cb\ 3> 1, we have: 

x (cb) = e~ (-cb - cb!) 1 -— + c t e* + - * , ■ (C2) 

(l + e) J 4> 



In the evaluation of the inverse Laplace transformation (|66|) one can neglect two last terms in Eq. ( |C2[ ) if one considers 
the DF ft (y) at y close enough to 1. Substituting x (cb) into Eq.(B6h, rescaling integration variable as: 



—cb — cb\ — e Gl y « u = V u , 

a 

and taking into account a< 1, which leads to the following expression in the exponent in Eq. (|gg|) 

-6 + yx (cb) =61+ 1 + € Wu (l - u~~^) w 0! + M^wlnu , 
a 



one immediately arrives at the formulas ( |73| - 

To evaluate S (W), it is convenient to use the integration contour 9 (ulnu) = 0, or, introducing polar coordinates 
u = r exp (16): 

r = exp (-9 cot 9) . (C3) 

Substituting Eq.(|C3|) into Eq.@, we have: 

S( W ) = ^-V(9)e~ wv ^, V(6) = ±- e -*^ . (C4) 
./_„. 27r sin 9 
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The value of the integral may be estimated by the saddle-point method, looking for the maxima of the expression 
In V (9) — WV (9). The stationary point equation is: 



V'{6) 



W 



1 



V(6) 



0. 



As If > 1, the stationary point is 9 S = 0, and the asymptote of Eq.(C4) is: 

S (W) w {2ireW)~ 1/2 exp (-W/e) . 



(C5) 



(C6) 



On the other hand, if W <C 1, we have two stationary points ±9 S , V (±6* s ) = 1/W, 9 S = it — S, and for 8 <C 1 the 
following equation may be obtained: 



7T /7T 



1 



1 



the solution for which is: 

The asymptotics of S (W) turns out to be: 

S(W) 



n/8 



In (e/W) 
In In (e/W) 



In In (e/W) 



eW 



In (e/W) 



(C7) 
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l 
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^exp(-i^) 
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1.303 
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2.00 
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TABLE I. Threshold values ct for bond (BP) and continuous (CP) percolation and critical exponents evaluated by the 
RGMK method of this work compared to the exact (if known) or best possible numeric values. Values for comparison were 
extracted from reviews [H,p2|, for the conductivity exponent fi after Ref. Ej|. 
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